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ABSTRACT 

We present the discovery of a strongly phase-variable absorption feature in the X-ray spectrum 
of the nearby, thermally-emitting, isolated neutron star RX J0720.4—3125. The absorption line was 
detected performing detailed phase-resolved spectroscopy in 20 XMM-Newton observations, covering 
the period May 2000 - September 2012. The feature has an energy of ^750 eV, an equivalent width of 
^ 30 eV, and it is significantly detected for only ^ 20% of the pulsar rotation. The absorption feature 
appears to be stable over the timespan covered by the observations. Given its strong dependence on 
the pulsar rotational phase and its narrow width, a plausible interpretation is in terms of resonant 
proton cyclotron absorption/scattering in a confined magnetic structure very close to the neutron star 
surface. The inferred field in such a magnetic loop is Hioop ^ 2 x 10^“^ G, a factor of ~ 7 higher than 
the surface dipolar magnetic field. 

Subject headings: X-rays: stars — stars: neutron — stars: individual (RX J0720.4—3125) 


1. INTRODUCTION 

RX J0720.4—3125 belongs to a group of seven 
thermally-emitting, radio-quiet, and nearby (;< 500 pc) 
isolated neutron stars, originally discovered in ROSAT 
all-sky survey data, and often referred to as the X-ray 
Dim Isolated Neutron Stars (XDINSs). They are among 
the closest known neutron stars, and are characterized by 
X-ray luminosities Lx ~ 10^°“^^ergs“^, long spin peri¬ 
ods (P ^ 3-11 s) and inferred surface dipolar magnetic 
fields Pdip ~ 10^^ G, partially overlapping those of the 
magnetars (see van Kerkwijk & Kaplan 2007 and Turolla 
2009 for reviews). XDINSs have estimated ages of a few 
10^ yr, derived from cooling curves and kinematics, while 
the characteristic ages are somewhat longer (~10®yr). 
The X-ray spectra show thermal emission with inferred 
temperatures kT « 50-100 eV and lack the non-thermal, 
power-law component often observed in other isolated 
neutron stars. Although a blackbody model provides an 
overall good description of their X-ray spectra, broad ab¬ 
sorption features have been found in six of the XDINSs 
with XMM-Newton observations. Their origin is unclear: 
they can be produced either by proton cyclotron reso¬ 
nances/atomic transitions in a magnetized atmosphere, 
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or by an inhomogeneous surface temperature distribu¬ 
tion, as recently suggested by Viganb et al. (2014). 

RX J0720.4—3125, the second brightest member of the 
class, was discovered by Haberl et al. (1997) as an iso¬ 
lated, pulsating neutron star with an 8.39 s spin period. 
The X-ray spectrum is best modelled by a blackbody 
{kT ^ 85 eV) plus a broad (~ 70 eV wide) absorption 
feature centered at ~ 270 eV (Haberl et al. 2004). A fur¬ 
ther, narrow (~ 5 eV wide) line was identified in XMM- 
Newton Reflection Grating Spectrometer (RGS) data at 
^ 570eV (Hohle et al. 2012a and references therein), 
possibly due to Oxygen circumstellar/interstellar absorp¬ 
tion. Unlike other XDINSs, RX J0720.4—3125 exhibits 
long-term variations in its timing and spectral parame¬ 
ters. In the period 2001-2003 the total flux stayed almost 
constant, whereas the blackbody temperature increased 
from ~ 84 to over 94 eV, and also the pulse profile, the 
blackbody radius, and the equivalent width of the ab¬ 
sorption feature changed. In the following years this 
trend appeared to reverse, with a decrease of the surface 
temperature. Two main explanations for this behaviour 
were proposed: either free precession, which predicts a 
cyclic pattern (Haberl et al. 2006), or a glitch that oc¬ 
curred around MJD = 52866±73days (van Kerkwijk et 
al. 2007). According to the most recent analyses (Hohle 
et al. 2009, 2012b), the interpretation in terms of a sin¬ 
gle sudden event around the proposed glitch epoch seems 
favoured. 

In this Letter we reanalyze all the archival XMM- 
Newton observations of RX J0720.4—3125, performing a 
detailed phase-resolved spectral analysis. In Sectionj^we 
report on the data analysis, and present the discovery of 
a second, new, phase-dependent absorption feature in the 
X-ray spectrum, the possible origin of which is discussed 
in Section [3l 

2. XMM-NEWTON DATA ANALYSIS AND RESULTS 

RX J0720.4—3125 was observed 20 times by the XMM- 
Newton satellite using the European Photon Imaging 
Camera (EPIC). We used here only data obtained with 
the EPlC-pn camera (Striider et al. 2001) because less 
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TABLE 1 

Summary of the XMM-Newton/EPlC-PN observations of RX J0720.4—3125*^ 


Obs. ID 

Obs. Date 
YYYY-MM-DD 

Read-out mode / filter 

Live time 
(ks) 

Source net count rate 
(counts s“^) 

Pile-up fraction ratios 
r=0”-30” 

0124100101 

2000 May 13 

FF / thin 

42.8 

6.46(1) 

0.963(3) 

0132520301 

2000 Nov 21 

FF / medium 

22.7 

5.60(2) 

0.964(4) 

0156960201 

2002 Nov 06 

FF / thin 

25.6 

6.60(2) 

0.969(3) 

0156960401 

2002 Nov 08 

FF / thin 

27.1 

6.54(2) 

0.966(3) 

0158360201 

2003 May 02 

SW / thick 

51.0 

3.480(8) 

1.011(3) 

0161960201 

2003 Oct 27 

SW / thin 

12.6 

7.52(2) 

1.013(5) 

0164560501 

2004 May 22 

FF / thin 

32.0 

6.96(1) 

0.971(3) 

0300520201 

2005 Apr 28 

FF / thin 

38.1 

6.86(1) 

0.968(3) 

0300520301 

2005 Sep 22 

FF / thin 

39.1 

6.93(1) 

0.969(3) 

0311590101 

2005 Nov 12 

FF / thin 

33.5 

6.75(1) 

0.970(3) 

0400140301 

2006 May 22 

FF / thin 

17.6 

6.83(2) 

0.970(4) 

0400140401 

2006 Nov 05 

FF / thin 

17.6 

6.90(2) 

0.966(4) 

0502710201 

2007 May 05 

FF / thin 

17.4 

6.80(2) 

0.968(4) 

0502710301 

2007 Nov 17 

FF / thin 

20.1 

7.71(2) 

0.971(4) 

0554510101 

2009 Mar 21 

FF / thin 

16.7 

6.84(2) 

0.967(4) 

0601170301 

2009 Sep 22 

FF / thin 

15.0 

6.77(2) 

0.968(4) 

0650920101 

2011 Apr 11 

FF / thin 

17.6 

6.61(2) 

0.973(4) 

0670700201 

2011 May 02 

FF / thin 

23.6 

6.73(2) 

0.965(3) 

0670700301 

2011 Oct 01 

FF / thin 

22.2 

6.60(2) 

0.972(3) 

0690070201 

2012 Sep 18 

FF / thin 

22.3 

6.60(2) 

0.970(3) 


^ FF: full-frame (time resolution of 73 ms); SW: small window (time resolution of 6 ms). Live time refers to the duration of the 
observations after filtering for background flares (see text). Count rates refer to the spectra extracted within a circular region with 
PATTERN = 0. Errors on the count rates are quoted at the Icr confidence level. Pile-up fraction ratios were calculated for single 
events alone and in the 0.1-1.2 keV energy range using the SAS epatplot tool. 
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Fig. 1.— Top panel: pulse profile for the May 2nd, 2003 ob¬ 
servation (red) and for all observations merged together (black). 
Bottom panel: normalized energy versus phase image obtained by 
binning the EPIC source counts into 0.01 phase bins and 25 eV 
energy channels for the observation performed on May 2nd, 2003. 


affected by pile-up (see Table [T|. We processed the 
data using the Science Analysis software (SAS, version 
13.5.0), adopting the most recent calibration files avail¬ 
able, and we removed any particle flare via good-time- 
intervals. All arrival times were referred to the Solar Sys¬ 
tem barycenter (source coordinates RA = 07^20“24!961, 
Dec = -31°25'5(y.'21, J2000.0; Kaplan et al. 2003), and 
a rotational phase was assigned to the source counts of all 
the observations using the timing solution of Hohle et al. 
(2012b). For the source distance we assumed D = 286 pc 
(Tetzlaff et al. 2011). All errors are reported at the 90% 
confidence level. 


Spectral analysis was performed using the XSPEC anal¬ 
ysis package (version 12.8.2, Arnaud 1996), using a min¬ 
imum of 100 counts per bin for each spectrum, the 
statistics and a maximum oversampling of the spectral 
energy resolution of a factor of 3. The fits were restricted 
to energies between 0.3 and 1.2 keV to exclude the broad 
absorption line at ~270eV so as to reduce the number 
of degrees of freedom (dof). We extracted the source 
spectra from a circular region of radius 30" centered on 
the source point spread function (PSF) and the back¬ 
ground counts from a circle of the same size far from the 
source and on the same CCD. For Small Window and 
Full Frame observations the pile-up level, estimated by 
means of the SAS epatplot tool, was <1.5% and 4%, 
respectively (see Table . To mitigate the pile-up we 
restricted our spectral analysis to photons having flag 
= 0 and PATTERN = 0. 

We started our analysis from the May 2nd, 2003 ob¬ 
servation, this being the longest observation, and least 
affected by pile-up (the pn was operating in Small Win¬ 
dow mode). The phase-averaged spectrum is well fitted 
by an absorbed^ blackbody (reduced chi-square xi = 
1.17 for 130 dofjyielding the following parameters: col¬ 
umn density Ah = 1.9(3) x 10^°cm“^, blackbody tem¬ 
perature kT^B = 82.1±0.f3eV and blackbody radius Rbb 
= 6.3±0.6km. 

To study in detail the spectral variability with phase, 
we produced a phase-energy image by binning the EPIC- 
pn source counts into 100 rotational phase channels and 
25-eV-wide energy channels, and then normalizing to the 
phase-averaged energy spectrum. The image (Fig. 

We took into account the effects of interstellar absorption 
along the line of sight through the TBABS model with the vern cross 
sections (Verner et al. 1996) and the wilm abundances (Wilms et 
al. 2000). 
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Fig. 2.— From top to bottom. Panel a: spectrum of the 0.1-0.3 phase interval (black circles) fitted with an absorbed blackbody plus 
a Gaussian profile (gauss model); the spectrum in the phase range 0.5-0.7 together with the best fitting model (tbabs*bbodyrad) is also 
shown (red triangles). Both spectra are from the May 2nd, 2003 observation. Panel b: residuals with respect to these models. Panel c: 
residuals of the previous spectra after setting the line normalization to zero. Panel d: residuals of the spectrum relative to phase 0.1-0.3 
after merging the phase-resolved spectra of all observations listed in Tablet. 


shows a feature in the phase interval 0.1-0.3, which gives 
a strong hint for the presence of an absorption line at 
^750 eV in the corresponding spectrum. Above ^ 1 keV 
the source counts become background dominated. 

To better investigate the significance of the feature ap¬ 
pearing in Fig. 1, we performed a phase-resolved spectral 
analysis for the May 2nd, 2003 observation dividing the 
rotational cycle in five equal phase bins. For the spec¬ 
trum relative to phase 0.1-0.3 the addition of a Gaus¬ 
sian line in absorption to the blackbody continuum leads 
to a significant improvement in the shape of residuals, 
while a simple blackbody model gives an acceptable fit 
for all the other phase bins. Both the gabs and gauss 
models result in xl = 0-88 for 89 dof, corresponding to 
an F-test probability of 3.5x10“^. We multiplied this 
value by 5, considering the five phase-intervals as tri¬ 
als. In this way we obtained a significance of ~ 5a. 
The values for line energy and width, obtained from the 
two models, are compatible within the errors. The best- 
fit tbabs* (bbodyrad + gauss) model gave the follow¬ 
ing parameters: feTeB = 83.0±0.1 eV, i?BB = 6.3±0.3 
km, line energy E;i„e = 745^27 eV, width a = 42^33 eV 
and normalization of -9.2lg'2 x 10“® (the column density 
was frozen at the value obtained for the phase-averaged 
spectrum). The equivalent width of the feature is 28)^11 
eV. Table 2 summarizes the results of the phase-resolved 
spectroscopy for the May 2nd, 2003 observation using the 


tbabs* (bbodyrad + gauss) model. Fig. shows the 
phase-resolved spectrum relative to phase O.PO.3 (black) 
in comparison with that at phase 0.5-0.7 (red): the ab¬ 
sorption feature clearly shows up in the residuals. 

The same phase-averaged and phase-resolved analysis 
was then repeated for all the remaining XMM-Newton 
observations in search for the presence of the 745 eV ab¬ 
sorption line. Since observations in Full Frame mode 
were affected by pile-up, their spectra were extracted 
from an annular region extending from 15" to 30" from 
the nominal source position, this time including also 
events with PATTERN<4. The resulting values of the 
line equivalent width and the F-test probability as a func¬ 
tion of time are shown in Fig. (top panel), and clearly 
indicate that the presence of the feature is unconstrained 
in most of the observations, possibly because of their 
lower statistics. 

We then performed a simultaneous fit of the spectra 
in the 0.1-0.3 phase bin for all the observations in Full 
Frame mode with the Thin filter (18 spectra in total; 
see Table 1). We constrained A^h and the line energy to 
be the same (to minimize the free parameters), and the 
width of the feature was frozen at the value obtained 
for the May 2003 observation. The phase-dependent 
feature is detected at an energy of 785(13) eV for the 
tbabs* (bbodyrad+gauss) model (x^ = 0.96 for 556 
dof); normalizations and equivalent widths are consistent 
within the errors in all the spectra. The tbabs*bbodyrad 
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TABLE 2 

Pulse phase spectroscopy for May 2nd, 2003 observation. 


Parameter" 

0.1-0.3 

0.3-0.5 

0.5-0.7 

0.7-0.9 

0.9-1.1 

BB 

kTBB (eV) 

81.1(6) 

81.9(6) 

82.2(6) 

82.6(6) 

82.5(9) 

Kbb (km) 

6 4^^'® 
'^•^-0.5 

6 

'^•^-0.5 

6 

'^•-‘--0.5 


6 

'^•^-0.6 

Flux'’ 

0.86(1) 

0.82(1) 

0.86(1) 

0.95(1) 

0.95(1) 

Unabs. Flux'' 

1.07 

1.01 

1.06 

1.17 

1.17 

NHP 

7.2x10“2 

5.9x10-1 

6.9x10-1 

4.2x10-'! 

2.7x10-1 

xl 

1.22 

0.96 

0.92 

1.26 

1.09 

dof 

92 

90 

92 

97 

75 

BB-I-GAUSS 

kTBB (eV) 

83.0ti:? 

82.4(9) 

QO 9 + U. / 

82.7(8) 

82 

i?BB (km) 

6.0(6) 

e-oirs 

6 1+^-® 
'^•-‘^-0.5 

e-atas 

6.4(6) 

Eline (eV) 

7451^^ 

745 

745 

745 

745 

<T '’(eV) 


41.7 

41.7 

41.7 

41.7 

|Norm| 

l x 10-5 

<2.7x10-5 

<1.7x10-5 

<2.9x10-5 

<2.4x10-5 

Eq. Width (eV) 

28t®4 

6t® 

<8 

<11 

<13 

F-test (10“®) 

3.5x10-^ 

232 

1000 

750 

>1000 

NHP 

7.9x10-1 

6.0x10-1 

6.6x10-1 

3.6x10-2 

2.5x10-1 

xl 

0.88 

0.95 

0.93 

1.27 

1.11 

dof 

89 

89 

91 

96 

74 


^ Nn was frozen at the value obtained for the phase averaged spectra: A^h = 1.9xl0^®cm 
^ Fluxes are calculated in the 0.3-1.2 keV energy range and in units of (10~^^ ergs~^cm“^) 

Line energy and width were frozen at the value obtained for the phase interval 0.1 - 0.3: = 745 eV and cr = 42 eV. 

^ NHP is the Null Hypothesis Probability. 




Fig. 3.— Top panel: evolution of the F-test probability and line 
equivalent width with time for the spectra in the 0.1-0.3 phase 
range for all XMM-Newton observations (the two observations per¬ 
formed during November 2002 were merged being only two days 
apart). Bottom panel: F-test probability and equivalent width as 
a function of the spin phase for the spectrum of all merged obser¬ 
vations. 


model gives xt = 1-03 for 575 dof. The F-test probabil¬ 
ity for the inclusion of the Gaussian absorption line is 
4.2x10“®, which corresponds to a significance level of 
^ 4.5cr (taking into account the number of trials). 

To further test the presence of this feature in data 
taken before and after May 2nd, 2003, we also combined 
all the spectra except that of May 2nd, 2003. We used 
the SAS epicspeccombine tool to merge the spectra and 
the response matrices. For the combined spectrum in 
the phase range 0.1-0.3, the best-fit with a blackbody 
model yields xt — 2-04 for 123 dof; by including a Gaus¬ 
sian absorption line, xt decreases to 1.48 for 121 dof. 
The parameters for the phase-dependent line are: 

= 795^^4 eV, normalization of -8.3^f 7 and equivalent 
width of 27^5 eV for the gauss model (the line width 
was fixed to the value obtained from the May 2nd, 2003 
observation). The F-test probability between these two 
models is 1.3x10“® and the significance is equal to 6a 
(again considering 5 trials). 

As a further check, we built a combined spectrum 
comprising all observations. This spectrum, created by 
merging all the spectra extracted from a 15" - 30" an¬ 
nular region with PATTERN<4, has an exposure time 
of 406.8 ks and a much improved statistics. The 0.1-0.3 
phase-resolved spectrum shows an absorption feature at 
energy 787^4® eV with normalization of -7.9^4 ® x 10“® 
and equivalent width of 26^5 eV for the gauss model; 
also in this case we froze the line width at the value 
estimated for the longest observation. The F-test prob¬ 
ability for the inclusion of the absorption feature is now 
4.7x10“® and, taking into consideration the number of 
trials, the significance level is again ^ 6a. As shown 
in Fig. panel d of, the spectral feature is clearly de- 
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tectable in the fit residuals of the merged observations, 
and the line significance increases with respect to the 
May 2003 observation alone. 

In the spectra relative to the other phase intervals the 
line equivalent width is always <13 eV^nd the F-test 
probability is never < 3 x 10“^ (see Fig. M lower panel). 
We stress that the response of the pn camera has slightly 
changed (see, e.g., Sartore et al. 2012) over the twelve 
years spanned by the observations, hence mixing spectra 
obtained with different settings over such a long period 
certainly introduce large systematic errors which are dif¬ 
ficult to quantify. However, a spectral feature present 
only in a given phase interval cannot result from sys¬ 
tematic effects, which are definitely independent on the 
rotational phase. While the line is mostly evident in the 
May 2nd, 2003 observation, which we used to derive the 
best estimate of the line parameters, the considerations 
presented above strongly suggest that the feature was 
present also in the other observations, with properties 
which appear relatively stable over the 12 yr period cov¬ 
ered by XMM-Newton data. 

3. DISCUSSION 

A careful re-analysis of all available XMM-Newton ob¬ 
servations of the isolated neutron star RX J0720.4—3125 
revealed the presence of a phase-dependent absorption 
feature at ^ 750 eV, which is mostly evident during the 
decline of the pulse profile (phases 0.1-0.3). The fea¬ 
ture was first signihcantly detected (at ^ 5cr confidence 
level) in the May 2nd, 2003 XMM-Newton observation, 
which was the longest one and was less affected by pile- 
up. Phase-dependent spectral analysis of the remaining 
datasets, which cover the period May 2000-September 
2012, indicates that the feature is likely present at all 
epochs and its properties are consistent with being con¬ 
stant over the timespan covered by available observations 
(^12 yrs). This suggests that the feature may be long- 
lasting and not associated to the timing anomalies (possi¬ 
bly due to a glitch) which occurred in 2003 (van Kerkwijk 
et al. 2007; Hohle et al. 2012b). 

Interestingly, an absorption feature with somewhat 
similar properties has been discovered in the spectrum 
of the “low-field” magnetar SGR 0418-1-5729 (Rdip 
6 X 10^^ G; Tiengo et al. 2013). In both cases the feature 
is present only in a quite narrow phase interval (^ 1/5 
of the phase cycle), but in SGR0418-1-5729 the line en¬ 
ergy is higher (> 2 keV) and strongly variable in phase 
(by a factor ~ 5 in one tenth of the phase cycle). No 
clear variation of the line energy with phase is seen in 
RX J0720.4—3125, although this might simply reflect the 
fact that the source counts become background domi¬ 
nated at energies just above that of the feature. Tiengo et 
al. (2013) interpreted the feature in SGR 0418-1-5729 as 
due to proton cyclotron resonant scattering in a baryon- 
loaded magnetic loop close to the star surface. The B- 
field variation along the loop explains the strong depen¬ 
dence of the line energy with phase, since photons emit¬ 
ted from the surface intersect the loop at different posi¬ 
tions. 

The similarities between the two features are sug¬ 
gestive that the same mechanism may be at work 
in RX J0720.4—3125. If the feature is a (broad¬ 
ened) proton cyclotron line, its energy is Aiine = 
0.63(i?ioop/10^^ G)/(l-|-z) keV, where 1/(1-|-z) ~ 0.7 is 


the gravitational redshift at the star surface. The implied 
magnetic held in the loop is Bioop ^ 1.8 x 10^^ G, a factor 
of ~ 7 larger than the the equatorial value of the dipole 
magnetic held at the surface, Rdip ~ 2.5 x 10^^ G. A pre¬ 
liminary calculation using the same geometrical model 
as in Tiengo et al. (2013) shows that the energy of 
the line and its appearance only for ^ 20% of the star 
period can indeed be reproduced. It should be noted, 
however, that the exact geometry of the magnetic struc¬ 
ture is likely different in the two sources. In the case 
of SGR 0418-1-5729, in fact. X-ray photons come from a 
quite small, hot spot with size ~ 1 km, while the emit¬ 
ting region of RX J0720.4—3125, albeit not uniform, is 
larger {Rbb ~ 6 km). 

A further possibility is that the line results from atomic 
absorption in the neutron star (magnetized) atmosphere. 
The strong dependence of the feature on the star rota¬ 
tional phase may be due to local inhomogeneities in the 
star surface magnetic field (as in the cyclotron scenario 
above) and/or to a clumpy structure of the atmosphere 
itself. The possibility that primary photons come from a 
small hot spot appears, on the other hand, unlikely, given 
the relatively large radiation radius and the broad shape 
of the pulse profile. If the feature arises from bound- 
bound (or bound-free) transitions, the atmosphere com¬ 
position is most probably He or mid-Z elements (see, 
e.g., Pavlov & Bezchastnov 2005; Mori & Ho 2007). The 
ionization energy for H is, in fact, below ~ 0.7 keV for 
fields < 10^"^ G (above this value the feature may be 
washed out by vacuum resonance mode conversion, e.g. 
Van Kerkwijk & Kaplan 2007 and references therein). 

The presence of proton cyclotron lines produced by 
resonant scattering/absorption in confined, high-R struc¬ 
tures close to the star surface, where the magnetic field 
is a factor > 10 higher than the dipole, would be sup¬ 
portive of a picture in which the magnetic field of (highly 
magnetized) neutron stars is complex, with substantial 
deviations from a pure dipole on the small scales. In 
recent years, a better (although far from conclusive) un¬ 
derstanding of how stellar magnetic fields are generated 
and amplihed indicates that strong, non-dipolar field 
components are ubiquitous, from massive stars (Braith- 
waite & Spruit 2005) to proto-neutron stars (Ober- 
gaulingiler, Janka & Aloy 2014). It is now well estab¬ 
lished that internal toroidal field components, multipolar 
surface structures, as well as very localized B-field bun¬ 
dles must be present in many neutron stars, especially in 
highly-magnetized ones. The recent discovery of sources 
showing strong magnetic-powered flares and outbursts, 
but with low dipolar magnetic fields, < 10^^ G, further 
strengthens this idea (Rea et al. 2010, 2012, 2014; Tur- 
olla & Esposito 2013). 

Finally, we stress that the absorption feature reported 
here for the first time is unlikely related to the broad 
absorption feature at 270eV (Haberl et al. 2004, 2006). 
Even if the two may have the same physical origin, i.e. 
a proton cyclotron resonance, the latter implies a lower 
magnetic field, ~ 5 x 10^^ G, comparable to the spin- 
down estimate, Rdip ~ 2.4 x 10^^ G, indicating that 
in this hypothesis it is related to the large-scale field. 
Actually broad absorption features have been detected in 
all the XDINSs, with the exception of RX J1856.5-3754 
(e.g. Turolla 2009), but the strong phase variability we 
observe for the 0.75 keV line of RX J0720.4—3125 sets it 
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apart. 

This is the first time evidence of a complex magnetic 
field structure is observed in an XDINS, although this 
result is not unexpected, if, as suggested by many in¬ 
vestigations, XDINSs are aged magnetars (Vigano et al. 
2013; see also Turolla 2009). 
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